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Abstract 

We study a system of nonlinear partial differential equations resulting 
from the traditional modelling of oil engineering within the framework 
of the mechanics of a continuous medium. Recent results on the problem 
provide existence, uniqueness and regularity of the optimal solution. Here 
we obtain the first necessary optimality conditions. 
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1 Introduction 

We are interested in the optimal control of the dead oil isotherm problem: 

'd t u - Atp{u) = div (g{u)Vp) in Q T = fl X (0, T) , 
d tP - div (d(u)Vp) = / in Q T = n x (0, T) , 

u \an = °> u \t=o = u o, 
Ulan = °. P\t=o=Po, 

where ft is an open bounded domain in R 2 with a sufficiently smooth bound- 
ary. Equations ([T]) serve as a model for an incompressible biphasic flow in a 
porous medium, with applications to the industry of exploitation of hydrocar- 
bons. To understand the optimal control problem we will consider here, some 
words about the recovery of hydrocarbons are in order. At the time of the first 
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run of a layer, the flow of the crude oil towards the surface is due to the energy 
stored in the gases under pressure in the natural hydraulic system. To mitigate 
the consecutive decline of production and the decomposition of the site, water 
injections are carried out, well before the normal exhaustion of the layer. The 
water is injected through wells with high pressure, by pumps specially drilled 
to this end. The pumps allow the displacement of the crude oil towards the 
wells of production. The wells must be judiciously distributed, which gives rise 
to a difficult problem of optimal control: how to choose the best installation 
sites of the production wells? This is precisely the question we address in this 
work: our main goal is to present a method to carry out the optimal control of 
(fT]) with respect to all the important parameters that intervene in the process. 
More precisely, wc seek necessary conditions for the admissible parameters u. p 
and / to minimize the functional 

J(u,pJ) = \ \\u-u\\1 Qt + \\\p-p\\1 Qt + ^ \\f\\Z, QT +Y W^Iqt ( 2 ) 

where qo > 1 and (3i > and /?2 > are two coefficients of penalization. The 
first two terms in (f5]) make possible to minimize the difference between the 
reduced saturation of oil u, the global pressure p and the given data U and P. 

Existence and uniqueness to the system (fl]) , for the case when the term dtp is 
missing but for more general boundary conditions, is established in pQ. In [2] we 
obtain conditions which provide existence and regularity of the optimal solutions 
to the problem of minimizing ([2|) subject to ([1]). Here we are interested to obtain 
necessary optimality conditions which permit to find the solutions predicted by 
the results in [2]. This is, to the best of our knowledge, an important open 
question. 

Several techniques for deriving optimality conditions are available in the 
literature of optimal control systems governed by partial differential equations 
[21 HI [5]. We obtain the optimality conditions by making use of a Lagrangian ap- 
proach recently used with success by O. Bodart, A. V. Boureau and R. Touzani 
for an optimal control problem of the induction heating [6], and by H.-C. Lee 
and T. Shilkin for the thermistor problem [7]. 

2 Preliminaries 

Our main objective is to obtain necessary conditions for a triple (u, p, f) to min- 
imize ([2]) among all the functions (u,p, /) verifying |T]). The intended necessary 
optimality conditions are proved in $3] under adequate hypotheses on the data 
of the problem, which provide regularity of the optimal solution. 

2.1 Notation and Functional Spaces 

In the sequel we assume that tp, g and d are real valued C 1 -functions satisfying: 
(HI) < ci < d(r), <p{r) < c 2 ; \d'(r)\, \tf/(r)\, \<p"(r)\ < c 3 Vr e R. 
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(H2) uo, po G C 2 (n), C/,P6 L 2 (Q T ), where u , p : ft -» K, U, P : Qt R, 
and u | afl = Polao = 0. 

We consider the following spaces: 

W*'°(Qt) := i p (0,T, w£(fi)) = { u g £ p (Qt), Vu G £ P (Q T )} , 

endowed with the norm \\u\\ w i,o {Qt) = \\u\\ p Qt + ||Vw|| P)Qt ; 

W^(Qt) ■= {u G W^°(Q T ), S7 2 u,d t u G £ p (Q T )} , 
with the norm \\u\\ w 2,i {Qt) = IMI^i.o^ + ||V 2 w|| p Qt + \\d t u\\ p QT ; 

V:={u£ Wl'°{Q T ), d t u G L 2 (0, T, W^fi)) } ; 

^ := {« e ^(Qt), ^ T=9 nx(o,T) = o} i 
T := {/ G L 2 "{Q T ), d t f e L 2 (Q T )} ; 

o 2-i 

H:=L 2 i(Q T )x W 2q 9 (fi). 
2.2 Coercive Estimate 

The following lemma provides a coercive estimate to linear parabolic systems 
that is useful for our purposes. 

Lemma 2.1 ( 8 ). Let fl be a bounded domain with a C 2 -boundary, and assume 
that 

A ijU £ C (Q T ) , b ijk G L r {Q T ) , cy G £ 5 (Qt) 
wzt/i r > n + 2 and ^4 satisfying the strong ellipticity condition: 

3 7o > : A ijk i(u)B i:i Bki > 7 o|B| 2 VB G M nXn , 
AijkiAkHj = AjiM = Aijik . 

Then, for any s G (l,f), s ^ %, and for arbitrary functions f £ L s (Qt,R™) 

and uo GT^ S (f2, K ) i/iere exists a unique solution u £ W 2,1 (Qt, R ) o/ i/ie 
proWem 

9 t Uj - A ijk i{x)u kj i + b ijk (x)u jk + Cij(x)uj = fi(x) , 
w| t=0 = u , u| ao = 0. 

Moreover, the estimate 

W U Ww^(Q T ) < c [\\f\\.,Q T + II M °H^-| (0) ) 

/io/rfs /or some constant c depending only on n, Q, T , 70 and the norms of the 
coefficients. 
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2.3 Existence of Optimal Solution 

The following existence theorem is proved in [5] using a technical lemma found 
in [5], Young's inequality and Aubin's Lemma, together with the theorem of 
Lebesgue and some compacity arguments of J. L. Lions [3]. The conclusion 
follows from the fact that J is lower semicontinuous with respect to the weak 
convergence. 

Theorem 2.2 ( 2J). Under the hypotheses (HI) and (H2) there exists a q > 
1, depending on the data of the problem, such that the problem of minimizing 
^ subject to ^ has an optimal solution (u,p,f) satisfying u £ W^'^(Qt), 
p G C([0,T];L 2 (tt)) nW^ q °(Q T ), d t p G L 2 (0,T, W^Sl)), f G L 2 «°(Qt), 
d t feL 2 (Q T ). 

2.4 Regularity of Solutions 

Regularity of the solutions given by Theorem 12.21 is also proved in [5j . Theo- 
rem [53] is obtained using Young's and Holder's inequalities, Gronwall Lemma, 
De Giorgi-Nash-Ladyzhenskaya-Uraltseva theorem, an estimate from [TU], and 
some technical lemmas found in 

Theorem 2.3 ( 2 ). Let (u,p, f) be an optimal solution to the problem of mini- 
mizing ^ subject to ([1]). Suppose that (HI) and (H2) are satisfied. Then, there 
exist a > such that the following regularity conditions hold: u G C a '% (Qt) , 
u,p£ W^°{Q T ), u,p£ W^(Q T ), d t u, d t p G i°° (0,T;L 2 (ft)) n W^°(Q T ), 

u e CHQ T ), u g W^ o (Qt), p g W*>l(Q T ). 
3 Main Results 

We define the following nonlinear operator corresponding to (p} : 

F :W xW xT — > H x H 
(u,p,f)—>F(u,p,f) = 

where 

, f \_( d t u- A(p(u) - div(g(u)Vp), j Q u-u 
*{u,p,j) y d t p-div(d(u)Vp)- f, 10 p-Po 

7o being the trace operator jqu = u\ t _ . Owing to the estimate 
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(see [TT]) and hypothesis (HI), we have 

tp"(u) \Vu\ 2 , g'(u)VuVp, d(u)VuVp G L^(Q T ) C L 2? (Qt) 
Thus, it follows that F is well defined. 
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3.1 Gateaux differentiability 

Theorem 3.1. In addition to the hypotheses (HI) and (H2), let us suppose that 
(H3) \<p"'\ < c. 

Then, the operator F is Gateaux differentiable and its derivative is given by 

SF{u,p,f)(e,w,h) = ^F(u + se,p + sw,f + sh) | s=0 = (SF U SF 2 ) 

_ I dte — div (</?'(w)Ve) — div (tp" (u)eVu) — div (g(u)Vw) — div (g'(u)e\/p) , 7oe 
Y dtw — div (d(u)\7w) — div (d'(u)eVp) — h, ^qw 

for all (e,w,h) £ W x W x T. 

Furthermore, for any optimal solution (u,p,f) of the problem of minimizing @ 
among all the functions (u,p, f) satisfying ([T]), the image of SF (u,p, /) is equal 
to H x H. 

To prove Theorem 13.11 we make use of the following lemma. 

Lemma 3.2. The operator SF(u,p, f) : W x W x T — > H x H is linear and 
bounded. 

Proof, of Lemma l3~2l We have for all (e, w, h) £ W x W x T 

S p F2(u,p, /)(e, w, h) = dtw — div (d(u)Vw) — div (d'(w)eVp) — h 
= d t w - d(u)Aw - d'(u)Vu.Vw - d'(u)eAp - d'{u)Ve.Vu - d\u)eVu.Vp - h , 

with S p F the Gateaux derivative of F with respect to p. Then, using hypothesis 
(HI), we obtain that 

\\5 p F 2 (u,p, f)(e,w,h)\\ 2q ,Q T < \\d t w\\ 2q ,Q T + ||Vw|| 2g ,Q T + c\\ Aw\\ 2q ,Q T 

+c||VM.Vw|| 2g ,Q T +c||eAp||2 9 ,Q r +c||Ve.Vu|| 2(? ,Q r +c||eVw.Vp||2 g ,Q T + ||/i||2 (Z ,Q T • 

(3) 

In what follows we consider the term ||eVit.Vp||2 g ,Q T . Similar arguments apply 
to the remaining terms of |3]). We have 

||eVw.Vp||2g,Q T < ||e||cx>,Q T ||Vw.Vp||2g,Q T 

< \\e\UQ T \\^u\\^_ QT \\Vp\U, QT 

< c\\u\\ w \\p\\ w \\e\\ w ■ 

Then, 

\\S p F 2 (u,p,f)(e,w,h)\\ 2q . QT < c{\\u\\ w , \\p\\ w , ||/|| T ) (Hellwr + \\w\\ w + \\h\\?) . 

(4) 
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On the other hand, 

6 u Fi(u,p, f)(e, w, h) = dte — div (ip (tt)Ve) — div (tp (u)eVtt) 
— div (g(u)Vw) — div (g (v,)eVp) 
= d t e - <p'(u)Ae - <p"(u)Vu.Ve 
- tp"(u)eAu - ip"(u)VeS7u - tp'" (u)e\Vu\ 2 
— g(u)Aw — g'(u)VM.Vw 

— g (u)eAp — </(u)Ve.Vp — g" (u)eVu.Vp , 

where S M F is the Gateaux derivative of F with respect to u. The same argument 
as above give that 

||*«Fi(u,p J /)(e,«; J fc)||a g ,g T < c(H|w, \\p\\w, ||/[[r) (||e|| w + Hk + ||fc|| T ) . 

(5) 

Consequently, by (j4]) and |5]) we can write 

\\6F(u,p,f)(e,w,h)\\ HxHx y <c(\\u\\ w , \\p\\ w , ||/|| T ) (||e|| w + \\w\\ w + ||ft|| r ) . 

□ 

Proof, of Theorem 13. II In order to show that the image of 6F(u,p, /) is equal 
to H x H, we need to prove that there exists a (w, e,h) G W x W x T such that 

c^e — div ((/j'(u)Ve) — div (p" (w)eVlf) — diu (j(m)Vw) — dii; (</'(tt)eVp) = a , 
9jw — diw (d(u)Vw) — dif (d'(u)eVp) — h = (3 , 

eU = °> e \t=o= b > 



<•)<> 



= 0, w| t=0 =a. 



(6) 

for any (a, a) and (/3, b) G i?. Writing the system ((6]) for h = as 

<9 t e - <p'(u)Ae - 2p"(u)Vu.Ve - p"(u)eAu - p"'{u)e\Vu\ 2 , 
—g(u)Aw — j'(S)Vu.Vw — g'(u)eAp — </(u)Vp.Ve — g" (u)e\7u.\7p = a , 
9tto - d{u)Aw ~ d'(u)S7u.Vw - d'(u)eAp - d'(w)VM.Ve - d'{u)eVu.Vp = f3 , 

Han = Q > H t=0 = a > 

(7) 

it follows from the regularity of the optimal solution (Theorem I2.3|) that 

v'\u)Au,<p"'(u)\Vu\ 2 ,g\u)Ap,g'\u)Vu.Vp,d\u)Ap,d'(u)\7u.Vp G L 2qo{QT) , 
(p'\u)Vu,g'(u)Vu,g\u)Vp,d\u)Vu G L 4?0 (<2t) • 

By Lemma 12.11 there exists a unique solution of the system , hence there 
exists a (e,w,0) verifying ©. We conclude that the image of 5F is equal to 
H x H. □ 
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3.2 Necessary Optimality Condition 

We consider the cost functional J: WxWxT^M.^ and the Lagrangian L 
defined by 



£ ( U ,P, f,Pi,ei,a, b) — J {u,p, /) + ( F(u,p, f), 



Pi a 

ei, b 

where the bracket (■, ■) denote the duality between H and H' . 

Theorem 3.3. Under hypotheses (Hl)-(HS), if (u,p, f) is an optimal solu- 
tion to the problem of minimizing © subject to ([Tj), then there exist functions 
(eT,Pi) G Wj' (Qr) x W 2 ' (Qt) satisfying the following conditions: 

dtei + div (ip (u)Vei) — d'(u)Vp.VpT — 99 (u)VtZ.Vei — g (u)Vf>.VeT. =u — U , 

ei Ian = ' ei I*=t = ^ ' 

<9tf7T + div (d(u) VpT) + diu (<7(u) VeTf) =p — P , 

Pilan^ ' Prl i= T = > 

at' 1 



/3 2 + 2g o /3 1 |7|^- 2 7^pr '■ / ' 



t=o 



at 



0. 



(8) 



t=T 



Proof. Let (it,p, /) be an optimal solution to the problem of minimizing (|2|) sub- 
ject to ([T]). It is well known (cf. e.g. [12]) that there exist Lagrange multipliers 
((pi, a), (el, b)) e H' x H' verifying 

8^ p j)£{u,p, /,pT,eT,a,6) (e,w,h) = V(e,w,h) eff xlfxT, 

with 8r UjV j\C the Gateaux derivative of £ with respect to (u,p,f). We then 
obtain 

((5 - f7)e + (p - P) w + 2<7o/3i |7| 290 " 2 /^ + fod~fd t h) dxdt 
(c\e — div ((p (u)Ve) — dif (y"(u)eVu) — dw (g(u)'Vw) — div (g' (u)eVp))eT dxdt 

(dtw — div (d(u)\7w) — div (d (u)eVp") — h)p\dxdt 

-<7 e,a) + -(7010,6) = V(e,tu,h) e W x W x T. 
This last system is equivalent to the following one: 

((w — C/)e — dw (d (it)eVj>) pi + c\e el — div (<p (u)Ve) el 

— dra (ip"(«)eVti) el — dit> (g'(u)eVp) eT) dxdi 

+ / (Q? — -P) 10 + dtwpT — div (d(u)\7w) pT — div (g(u)\7w) el) dxdt 



Qt 



+ 



(2q p 1 \f\ 2 ^- 2 fh + fcdtfdth - flh) dxdt 



+ (7 e, a) + (70W, b) = V(e, to, h) £ W x W x T. 
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In others words, we have 

((u -U) + d'(u)Vp.S7pT - d{el - div (<p'(u)Vei) + f" (u)Vu.Vel + g'(u)Vp.Vei) e dxdt 

+ / {{p — P) + d t pi — div (d(u)Vpi) — div (g(u)Vei)) w dxdt 
Jq t 

+ f (2q p 1 \f\ 2qo - 2 Jh + fcdtjdth - pTh) dxdt 
+ (7oe,a) + (7010, b) = QM{e,w,h) G W x W x T. 

(9) 

Consider now the system 

d t e x + div (yj'(u)Vei) - d'(u)Vp.Vpi - ^"(w)Vu.Ve! - g'(u)Vp.Vei =u-U , 
dtpi + div (d(u)Vpi) + efo (<?(u)Vei) =p — P, 
e ilan = Pi I an = °> e ilt=r = Pilt=T = °- 

(10) 

It follows by Lemma [2TTI that IjlOp has a unique solution (ei,pi) G W 2 ' (Qt) x 
^2' 1 (Qt')- Since the problem of finding (e,w) £ W x W satisfying 



dt& — div (^'(ti)Ve) — div (jp (u)eVu) — div (j(u)Vw) — diw (g'(u)eVp) = sign{e\ — e±) 
dtw — div (d(u)Viu) — (<i (u)eVp) = sign{p\ — pi) 
7 e = 7 w = 

(11) 

is uniquely solvable on x W 2 g by Lemma 12.11 choosing h = in Q , 

multiplying (|10p by (e, to), integrating by parts, and making the difference with 
©, we obtain 

(c\e — div {ip'{u)Ve) — div (tp" (u)eVu) — div (g(u)\7w) — div (g 1 (u)eVp)) (e\ ~ el) dxdt 
+ / (dtw — div (d(u)\/w) — div (d'(w)eVp)) (pi — pi) dxdt 

Jq t 

+ (7oe, 7oeI -a) + (70 to, 7 pT - 6) = V(e, w) G W x W. 

(12) 

Choosing (e, to) in (fT2)) as the solution of the system (jTTJ) , we have 

sign(ei — ei)(ei —eT) dxdt + / sign(pi — pT){pi — Pi) dxdt = . 

JQt 

It follows that ei = el and pi = pT. Coming back to (|12p . we obtain 7oeT = a 
and 70P1 = 5. On the other hand, choosing (e, w) = (0, 0) in (|9|) it follows (|8|), 
which conclude the proof of Theorem 13.31 □ 
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